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Abstract 

The conservation of the isovector axial current in lattice QCD with mass- 
less Wilson quarks is studied to one-loop order of perturbation theory. Fol- 
lowing a strategy described in a previous publication, the 0(a) counterterm 
required for on-shell improvement of the current is computed. We also confirm 
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improved action. 
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1. Introduction 



This paper is the second in a series of publications [1-3] on chiral symme- 
try and 0(a) improvement in lattice QCD with Wilson quarks. The underlying 
theoretical framework has been explained in detail in ref. [1]. We now report 
on a study of the conservation of the isovector axial current to one-loop order 
of perturbation theory. Our aim is to verify that the 0(a) improvement of the 
lattice action and the current works out in the expected way. In particular, we 
shall show that the coefficients of the associated 0(a) counterterms, c sw and 
ca, can be computed by requiring the lattice corrections to the PCAC relation 
to be of order a 2 (as suggested in ref. [1]). 

The coefficient c sw of the SW term [4] in the improved action has been 
calculated to one-loop order many years ago by Wohlert [5]. His calculation 
has partly been repeated by Naik [7], but it is only now that a complete and 
independent check of Wohlert 's result is obtained. The other coefficient, ca, 
vanishes at tree-level and a one-loop calculation is hence required for a first 
estimate of the magnitude of the associated 0(a) counterterm. 

Throughout this paper we assume that the reader is familiar with ref. [1]. 
The notations introduced there are taken over completely without further no- 
tice. Equations in ref. [1] are referred to by prefixing a Roman "I" to the 
equation number. 

The correlation functions involving the axial current that we shall consider 
are constructed from the Schrodinger functional. We first need to discuss the 
evaluation of the quark functional integral in this framework, which is not 
totally trivial because the boundary conditions are rather special (sect. 2). In 
sect. 3 the PCAC relation and the approach to the continuum limit are studied 
in the free quark theory. The perturbation expansion of the relevant correlation 
functions in the full theory is derived in sect. 4. We are then in a position to 
compute the one-loop coefficients c A *' and ciw (sects. 5,6). In the last section 
we show that the lattice effects that remain after improvement are rather small 
in general. 

We are indebted to our colleagues Andrea Galli, Stefan Sint, Rainer Som- 
mer and Ulli Wolff for their help at various stages of this work. We also thank 
the computer centers at CERN and EPFL for allowing us to use their SP2 and 
T3D machines. 



1 



2. Quark propagation 



In this section we consider the 0(a) improved Schrodinger functional de- 
fined in sects. 4 and 5 of ref. [1] and work out the integral over the quark and 
anti-quark fields in some detail. 

2.1 Lattice Dirac operator 

For all times xq in the range < xq < T we have 

S 4^2L = (D + 5D + m )^(x), (2.1) 
Sip(x) 

where D denotes the Wilson-Dirac operator (1.2.3). The 0(a) correction SD 
is a sum of a volume and a boundary term, 

SD = 5D V + 5D h . (2.2) 

Explicitly the volume term is given by 

SD v tjj(x) = c sw jaa llu F llu (x)il)(x), (2.3) 

while for the boundary term one obtains 

SD h ^(x) = (c t - l)^{^ ,a bl>(x) ~ U(x - a6, O^P+^x - a6)] 

+ 8 x<uT -a [${x) - U(x, 0)P_ij(x + a6)] }. (2.4) 

Note that in this expression there is no reference to the "undefined" compo- 
nents of the quark field at the boundaries of the lattice. D + SD is hence a 
well-defined operator acting on quark fields satisfying Schrodinger functional 
boundary conditions. 

2.2 Quark propagator and classical solutions 

In the presence of a given gauge field U(x, fj,), the quark propagator S(x,y) is 
defined through 

(D + 5D + m )S(x,y) = a- 4 5 xy , < x < T, (2.5) 
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and the boundary conditions 

P+S(x,y)\ Xo=0 = P_S{x,y)\ Xo=T = 0. (2.6) 

We assume that the gauge field and the quark mass are such that the solution 
of eq. (2.5) is unique. As in infinite volume one can show that 

S(x,yy = j 5 S(y,x)j 5 . (2.7) 

An analytic expression for the propagator in the free quark theory will be given 
in sect. 3. 

The solution of the Dirac equation 

(D + 8D + mo)i/) cl (x) = 0, < x < T, (2.8) 
with boundary values 

P + M*)\x =o = /f(x), P_VciML=T = /(x), (2.9) 
plays an important role in the following. For < xq < T it is given by 

4> c \(x) = a 3 ^2c t {S(x 1 y)U(y - a6 1 0)~ 1 p(y)\ yo=a + 
y 

S(x,y)U(y,0)p\y)\ yo=T -a}. (2-10) 

It should be emphasized that this representation is not valid if xq = or 
x = T. We shall also make use of the solution of the adjoint equation 

^ cl (a;)(S t + 5S t + mo) = 0, < x < T, (2.11) 
with boundary values 

^ cl (a;)P_U = o=/o(x), ^ cl (x)P + \ X0=T = p'(^). (2.12) 
The adjoint action of the Dirac operator in eq. (2.11) is defined through 

-flf2L = ^ (a;)( £t + ^ t + mo)) (2.13) 

and it is again possible to express the solution through the propagator in a 
way analogous to eq. (2.10). 
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2.3 Quark functional integral 

We are interested in evaluating the expectation value of polynomials O in the 
gauge field variables, the quark and anti-quark fields in the interior of the 
space-time volume and the boundary fields C(y), . . . , C'( z )- If we first integrate 
over the quark fields, the expectation value assumes the form 

(0) = ([0] F ) G , (2.14) 

where (. . .)q denotes the gauge field average with probability density propor- 
tional to 

det(D + 5D + m )exp{-S G [U] - 5S G)h [U]}. (2.15) 

The Dirac operator is here considered to be a linear mapping in the space of 
all quark fields with zero boundary values. 

In the following we deduce a set of simple rules that allow us to calculate 
the quark field average [0] F . We begin by introducing the generating functional 

Zf[P', p';p, P\% T,U] = J B[ip]D[ip] exp{-S F ,i mpr [t/,V>,V'] 

+ a 4 ^2l~ip(x)r](x) + f](x)'ip(x)]}. (2.16) 

X 

The improved quark action appearing in this formula is equal to the sum of 
all terms in the improved action (1.5.2) that depend on the quark fields. r](x) 
and f)(x), < xq < T, are source fields for the quark and anti-quark fields in 
the interior of the space-time volume. If we substitute 

* {x) ^-s£xj' (2 - 17) 

in the polynomial O, its quark field average is obtained through 

[0] F = \^-OZ F \ . (2.18) 

As will be shown below the generating functional is an exponential of a qua- 
dratic expression in the boundary values p' , . . . , p and the source fields 77,77. 
The differentiations in eq. (2.18) may hence be worked out by applying Wick's 
theorem, i.e. [0] F is expressed as a sum of Wick contractions. We are then 
left with the task to list all quark two-point functions. 
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To this end we insert the decomposition 



(2.19) 



in the functional integral (2.16). The quantum components of the quark fields, 
x(x) and x( x )i have vanishing boundary values. It is then easy to show that 

S^imprttA V', i>] = SF,impr[£A V'cl ) V'cl ] + S F)impT [U, X, X ] , (2.20) 

and for the generating functional we thus obtain 

lnZ F = In Z F \ pl = _ - v = ~ S'F.imprt^, V'cl) V'cl] 

+ a8 ^2v( x )S(x,y)i](y) + a 4 ^ [v(x)A\( x ) + ip cl (x)i](x)] . (2.21) 

x,y x 

The second term on the right hand side of this equation may be simplified 
by noting that (D + SD + mo)ip c i(x) vanishes for < xq < T. Taking the 
boundary conditions into account, some algebra then yields 

S^imprttA V'cl) V'cl] = 

a 3 J2\ \ a ~ c * t(x)7fc(V* + V fc )/o(x) + p'(x)y k (Vt + V fc )/o'(x) 



p{y)U{x - a6, O)^ci(^) L o=a + p'(x)U(x, 0)" Vci(z) L 0=T _ a 



(2.22) 



The basic two-point functions may now be calculated by differentiating with 
respect to the boundary values and the source fields. The complete list of 
non-zero contractions is 



[<ip(x)ip(y)] F = S(x,y), 



(2.23) 



[<%K(y)] F = c t S(x, y)U(y - a6, 0)"^+ \ yQ=a , 



[<%K'(y)] F = c t S(x, y)U(y, 0)P-\ yo=T _ a , 
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(2.24) 
(2.25) 



[C(x)?%)] F = c t P-f/^-aO, 0)5(^,^)1^^ (2.26) 

[C'(x)?(y)] F = c t P+C/^, O)" 1 ^, y) L =T _ a , (2.27) 

[C(x)C(y)] F = cl P_U(x - ad, 0)S(x, y)U(y - a6, 0)"^+ \ Xo=yo=a 

- ±c s P_ 7fc (v£ + V fc )a" 2 5 xy , (2.28) 

[C(x)C'(y)] F = cl P.U{x - a6, 0)S(x, y)U{y, Q)- l P-\ Xo=a , yo=T _ a , (2-29) 

[C'(x)C(y)] F = c t 2 P+f/(a ; ,0)- 1 5( a; ,y)f/(y-a6,0)- 1 P + | :co=T _ a ^ =a , (2.30) 

[C'(x)C'(y)] F = c t 2 P+f/(a ; ,0)- 1 5( a; ,y)f/(y,0)P_| :co=!/o=T _ a 

- ic s P +7fc (v£ + V fc )a" 2 5 xy . (2.31) 



We have here made use of eq. (2.10) and the analogous representation of ip c \{x) 
to express all contractions through the quark propagator. 

2.4 Example 

To illustrate the procedure we consider the correlation functions A( s o) an d 
fp(xo) introduced in subsect. 6.3 of ref. [1]. Applying Wick's theorem one 
deduces that 

AM = « 6 E i ( tr { [C(#)] F 7075 h%K(y)] F 7s}} G , (2.32) 

/p^o) = a 6 El(tr{ [C(z)#*)] F 75 ty(*)C(y)] F 75}> G , (2-33) 

where the contractions are given by eqs. (2. 24), (2. 26) and the trace is to be 
taken over the Dirac and colour indices only. Note that, as a consequence of 
the hermiticity property (2.7), we have 

75 [C(z)#r)] F 75 = {^WCW] F } f . (2.34) 

The correlation function fp(xo) is hence an average of a manifestly positive 
quantity. 
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3. Free quark theory 



We now proceed to study the propagation of free Wilson quarks with 
Schrodinger functional boundary conditions. The formalism developed above 
applies in this case too with the obvious modifications. We set c s = c t = 1 in 
this section and shall justify this choice in subsect. 3.3. 

3.1 Analytic expression for the quark propagator 

To introduce our notations we first discuss the plane wave solutions of Dirac's 
equation on lattices that are infinitely extended in the time direction. It suffices 
to consider the positive energy solutions 

^(x) = ue vpx , lmp >0. (3.1) 

The spatial components of the momentum p are in the range 

— ir/a < pk < it /a. (3-2) 

They must be integer multiples of 2ir/L to satisfy the boundary conditions in 
the space directions. The constant spinor u and the energy po are constrained 
by Dirac's equation, (D + mo)ip(x) = 0, which reduces to 

{ij fl pt + M(p + )}u = 0. (3.3) 

Here and below we use the abbreviations 

q^ = (I /a) sin (aq^) , (3.4) 

g M = (2/a) sin , (3.5) 

for any momentum q. The effective mass in eq. (3.3) is given by 

M{q) = too + \aq 2 , (3.6) 

and the momentum p^ = p^ + includes a shift by the angles 9^ (cf. sub- 

sect. 4.2 and appendix A of ref. [1]). 

Eq. (3.3) is similar to Dirac's equation in the continuum theory and its 
solution is obtained in essentially the same manner. We first multiply the 
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equation from the left with ij^p^ — M(p + ) and deduce that 

(p+) 2 + M(p+) 2 = 0. (3.7) 
Solving for the energy po one obtains 

Po = Po = ^(p + ) mod 2k /a, (3-8) 

where oj(q) is given by 

I 

a ( q 2 + (to + 7r«q 2 ) 2 ^ 2 



sinh 



2 UW 



2^1 + a(m + |aq 2 ^ 



(3.9) 



w(p + ) is well-defined and non-negative for too > 0. There are no other positive 
energy solutions of eq. (3.7). 

It remains to determine the spinor u. With po as above, eq. (3.3) has two 
linearly independent solutions, corresponding to spin up and spin down states. 
Explicit expressions could be given as in the continuum theory, but we shall 
not need them here. 

The propagator S(x,y) can now be constructed following standard proce- 
dures. We omit the details and state the result in the form 

S(x,y) = (D^ + m )G(x,y), 0<x o ,y o <T, (3.10) 

where G(x,y) is defined through 

G(x, y) = L~ 3 {-2^Po + A(p + )i?(p + )}- 1 e^ x -^ 
p 

X j(M(p+) - ipZ)e-^ p+ ^ Xo - y °\ + (M(p+) + ip+) e -"(p + )(2T-\ Xo -y \) 

- (M(p + ) + i 7o p f )e- w(p+)(a;o+!/o) - (M(p+) - i 7o p f )e- w(p+)(2T - a;o - !/o) }, 

(3.11) 

for all times xq (including and T). The amplitudes A and R occurring in 
this expression are given by 

A(q) = l + a(TO + iaq 2 ), (3.12) 



R(q) = M{q) jl - e" 2 ^ 7 } - iq jl + e" 2 ^ 7 } . (3.13) 

The sum runs over all momenta p in the range (3.2) and p + is on the mass- 
shell (3.8). Note that R(p + ) vanishes if p + = and too = 0. The propagator 
remains well-defined, however, and can be calculated by first assuming too to 
be positive and then taking the limit too — > 0. 

3.2 Axial current conservation 

As explained in sect. 6 of ref. [1] the PC AC relation in the improved theory 
is violated by cutoff effects of order a 2 . Without improvement the effects 
are of order a in general. The PCAC relation may hence be used as a test for 
improvement. Specifically it has been proposed to introduce an unrenormalized 
current quark mass to through eq. (1.6.13). Any variation of to as a function 
of the kinematical parameters T, L, etc., can be taken as a measure for the size 
of the lattice effects at the chosen lattice spacing. 

In the free quark theory the 0(a) counterterm 8A a ^(x) to the axial current 
is expected to vanish [6] and to is then given by 

m=\{dt + d )f A {x )/f ¥ {x ). (3.14) 

We are now in a position to compute this ratio analytically. Combining 
eqs. (2. 32), (2. 33) with (2. 24), (2. 26) and the exact expression for the free prop- 
agator, one obtains (after some algebra) 

/a(*o) = ^{2M + (/)M_(/)e-Mp + )T 

-M(p+) ^M-{p + )e- 2u ^ Xo +M + (p + )e- 2w ( p+ '( 2T - a;o )]}, (3.15) 

f P ( Xo ) = _||^{M_(p+)e- 2 -(P + )- - M + (p+)e- 2 -(P + M 2T --)}, (3.16) 

where N denotes the number of colours, p = and 

M±(p + ) = M{p + )±ipZ (3.17) 
[as before the momentum p + is assumed to be on the mass-shell (3.8)]. From 
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these formulae one immediately deduces that 

to = M(p + ) cosh[aw(p + )] (3.18) 

for all xq in the range a < xq < T — a. 

It may come as a surprise that to turns out to be independent of xq. Al- 
though this is the expected behaviour in the continuum limit, lattice effects in 
general imply that to is a non-trivial function of xq. If there were no lattice 
effects at all, to would also have to be independent of the kinematical parame- 
ters such as L and the angles 9k- This is evidently not the case, but from the 
above it follows that 

TO = TO p +0(a 2 ), (3.19) 

where to p denotes the pole mass (1.3.6). In other words, there is no dependence 
on the kinematical parameters at order a. 

3.3 Continuum limit and 0(a) improvement 

Since all correlation functions are obtained as products of the basic two-point 
functions (2.23)-(2.31), it suffices to study the approach to the continuum limit 
of the latter. Before taking a — > we must specify in which way the quark mass 
is to be treated. Our discussion above suggests that we should keep the pole 
mass TO p fixed instead of the bare mass too. The name "pole mass" derives 
from the observation that 

w(0) = \m p \. (3.20) 

For general momenta q we have 

^(q) = e(q)+0(a 2 ), e(q) = ^TO2+q2. (3.21) 

It is also straightforward to show that 

M(p + ) = to p + 0(a 2 ), (3.22) 

R(p+) = TO p {l - e - 2e (P + ' T } + e(p + ) {l + e - 2e (P + ' T } + 0(a 2 ), (3.23) 

for any momentum p + on the mass-shell (3.8). Using the explicit expression 
for the quark propagator, the continuum limit of the contractions (2.23)-(2.31) 
is now easily worked out. 
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We first consider the two-point functions of the boundary fields, 



a 3 £ e -*P(x- y) (C(x)C - (y)> = _ m p t |i£J {! _ 6 -Mp + )t} p +> (3 . 24) 

a 3^ e - 8P (x- y) ( C(x)C -' (y )) = -2.Po + ^|e-^ + ) T P_. (3.25) 

Since A(p + ) = 1 + arriQ + 0(a 2 ), it follows that these correlation functions are 
0(a) improved up to an overall normalization factor equal to 1 + arriQ. Note 
that the second term in eq. (2.28) cancels a momentum dependent contribution 
of order a. Without this term the correlation function (3.24) would not be 
improved. 

The mixed two-point functions, involving a boundary field and the bulk 
field ip(x) or ip(x), are also 0(a) improved. This is immediately clear from 

« 3 E e " 8p(x " y) <<%K(y)} = ^ry {(M(p + ) - tp + o - t7 k pt)e-^ x ° 

- (M(p+) + iffi - i lk pt)^ u( ^ )(2T ~ Xo) }P+- (3.26) 

No renormalization factor of the form 1 + O(amo) is found here, because the 
factors associated with the boundary and bulk fields are inverse to each other 
and thus cancel in the mixed two-point functions. An 0(a) normalization factor 
equal to (1 + amo) -1 in fact shows up in the correlation function (ip(x)ip(y)), 
as one may easily verify from eqs. (3. 10), (3. 11). 

We thus conclude that the continuum limit of the free quark Schrodinger 
functional is reached with a rate proportional to a 2 , provided the fields are 
renormalized by factors of the form 1 + b(gQ)am q and if we keep the pole mass 
nip fixed instead of the bare mass (cf. sects. 2,3 of ref. [1]). In particular, 
the PCAC relation holds up to corrections of order a 2 . If we expand the b- 
coefficients according to eq. (1.3.17), our calculations imply that 

b^=b^ = -\. (3.27) 

We can also easily take the continuum limit of /a^o) an d fp( x o) an d verify 
that these correlation functions are improved if we set 

&i 0) = 6< 0) = 1, (3.28) 
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thus confirming a result originally obtained in ref. [6]. 

As anticipated the 0(a) boundary counterterms in the improved quark ac- 
tion are not needed for 0(a) improvement at tree-level. One may in fact prove 
that they must be set to zero as otherwise one would end up with uncancelled 
0(a) terms in some of the correlation functions considered above. The same 
is true for the 0(a) correction to the axial current, i.e. the coefficient ca must 
vanish at go = 0. 



To determine the coefficients ca and c sw to one-loop order of perturbation 
theory we need to compute the correlation functions /a(#o) an d fp( x o) to 
order g% for various choices of the kinematical parameters. In this section we 
describe in some detail how this calculation is done for the case of vanishing 
boundary values C and C of the gauge field. The modifications required for 
non-zero boundary values will be discussed in sect. 6. 

With little additional effort we shall also be able to compute the one- 
loop coefficient cj 1 '. The associated 0(a) boundary counterterm is hence not 
dropped here, although we shall not need to include it when calculating the 
other coefficients (cf. eq. (1.5.25) and subsect. 6.1 of ref. [1]). 

4-1 Preliminaries 

To evaluate the expectation values (2.32) and (2.33) in perturbation theory, we 
first fix the gauge and then expand the action and the quark field contractions 
in powers of the coupling. As for the action we shall not need to go beyond the 
leading term. The expansion of the two-point contractions can be simplified 
by noting that 



4. Expansion of /a^o) an d fp( x o) to one- loop order 





So if we define the matrix 





and recall eq. (2.34), it follows that 




(4.3) 
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with Ta = —70 and Tp = 1. It is, therefore, sufficient to work out the pertur- 
bation expansion of ip c \(x) to order g$. 

4-2 Gauge fixing and the gluon propagator 

The calculation of the expectation value in eq. (4.3) involves an integration 
over all lattice gauge fields with the specified boundary values. Under a gauge 
transformation 

U(x,n) Q{x)U{x,ii)Q{x + ajl)- 1 , Q(x)eS\J(N), (4.4) 

the integrand and the boundary values are invariant, provided the gauge func- 
tion Q(x) is independent of x at the boundaries xq = and xq = T. 

To fix this gauge degeneracy we may follow the steps taken in sect. 6 
of ref. [9] for the analogous case of the Schrodinger functional in pure gauge 
theories. A subtle point to note is that the boundary conditions on the ghost 
fields c(x) and c(x) are not the same for all boundary values of the gauge field. 
For C = C = the boundary conditions are that the ghost fields should vanish 
at xq = T, while at xq = they may be non-zero but should be independent 
of x. The linear space of all fields with these boundary conditions [and taking 
values in the Lie algebra of SU(7V)] is denoted by C 

In the gauge fixed theory the link variables are parametrized through a 
vector field q^{x) according to 

U(x,fi) = exp{g aq l _ t (x)} . (4.5) 

Lattice vector fields q^{x) are defined for < xq < T (if fj, = 0) and < xq < T 
(otherwise). They take values in the Lie algebra of SU(7V) and satisfy the 
boundary conditions 

qk{x)\ Xo =o = qk{x)\ Xo =T = 0. (4.6) 

The space of all these fields is denoted by H. 

To be able to write the gauge fixing term in a compact form we introduce 
the operator 

d: C^Ti, (du) tl (x) = d tl u(x). (4.7) 

The adjoint d* of d maps any vector field q £ H onto an element of L such 
that 

(d*q,uj) = -(q,duj) for all u G C. (4.8) 
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With the obvious choice for the scalar products on H and L we have 

' d*q^(x) if < x < T, 

d*q(x) = | (a 2 /L 3 ) E y go (y) L=o if x = 0, (4.9) 
,0 if xq = T. 

Note that d*q(x) satisfies the boundary conditions required for fields in C 

The total action of the gauge fixed theory, including the gauge fixing term 
and the ghost field action, is now given by 

S tot [q, c, c] = S G [U] + 5S G)h [U] - tr {ln{D + 5D + m )} 

+ ±\ (d*q,d*q)-(c,d*5 c q), (4.10) 

where Ao denotes the bare gauge parameter and S c q the first-order variation of 
q under the gauge transformation generated by the ghost field c. We shall not 
need the explicit form for S c q here, but note for completeness that 

S c qf, = d^c + g [qv,c+ \ad^c] + 0(gl) (4.11) 

(no sum over fj, is implied in this formula) . In the following expectation values 
in the gauge fixed theory, with action (4.10) and a priori measure D[f/]D[c]D[c], 
are denoted by (. . .) G . 

To leading order in the gauge coupling we have 

Stot[q, c, c] = constant + \ (q, A t q) + (c, A c) + O(g ), ( 4 -12) 

where Ao = —d*d and 

Ai = A; - \ dd*, (4.13) 

A^q^x) = ^ {-dtd v q ll {x) + dtd^q^x)} . (4.14) 

The action of Ai simplifies considerably in the Feynman gauge Ao = 1. Ex- 
plicitly, for the spatial components one finds 

A 1 q k (x) = -dtd fl q k (x), 0<x o <T, (4.15) 
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while for the time component we have 



Aigo(z) = < 



if < xq < T — a, 
dtdkqo(x) + a~ l dtq (x) ifx =T-a, 
dtdkqo(x) ~ a~ l d q (x) 



(4.16) 



I +(«/£ 3 )E y goG/)| 



2/0 = 



if 



0. 



It is then straightforward to calculate the gluon propagator 

(q a , {x)q h M)G = S ab D^( Xl y) + (% 2 ). 
The result is given in appendix A. 



(4.17) 



4-3 Expansion of ip c \(x) in powers of go 

Through eq. (4.5) the Wilson-Dirac operator becomes a function of the bare 
coupling go- The lowest order term in the expansion 



D = J2g£D 



(k) 



(4.18) 



fc=0 



coincides with the free quark Dirac operator studied in sect. 3. For k > 1 we 
have 



D (k ^{x) 



,k-l 



k\ 



11=0 



+ i(l + l^)\~ l [-q^{x - afi)] k i;(x - afi)Y (4.19) 

The 0(a) counterterm SD can be expanded similarly. In this case the depen- 
dence on </o is also through the coefficients c sw and c t . Note that SD is of order 
go, because the gluon field tensor vanishes proportionally to the coupling and 
because c t — 1 is of order g$ (cf. sect. 3). 

The classical solution ip c \(x) may now be obtained in the form of a series 



(4.20) 



k=0 
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where the leading term, ip ci ( x )i satisfies the free Wilson-Dirac equation with 
boundary values p(x.) and /o'(x) (cf. subsect. 2.2). For < xq < T the solution 
is given by 

(x) = a 3 J2 {S {0) (*> V)P(V) L=- + S<°> (a, y)p'(y) |, 0=T -a } , (4.21) 
y 

with S*' ' (a:, y) being the free quark propagator discussed in sect. 3. The higher 
order terms have zero boundary values and may be recursively constructed by 
solving 

k-l 

(D^ +m )^ ) (x) = -J2(D {k ~ j) + SD^-^)^ ] (x) (4.22) 

i=o 

for < xq < T. In particular, to first and second order we have 

^(x) = -a^s (0) ( X ,y)D (l) ^(y), 

y 
y 

+ a s J2S {0) ( X ,y)D {1) S^(y,z)DW4 ( l\z), 

y,z 

where the abbreviation = + 8D^ has been used. 

4-4 Expansion of fx(xo) and Feynman diagrams 

The perturbation expansion of the matrix H(x) defined in subsect. 4.1 may be 
obtained straightforwardly by differentiating the series (4.20) with respect to 
the boundary value p(y). To lowest order we have 

and the higher order terms, (x) and H^(x), are given by eqs. (4.23) and 
(4.24) with ^(x) replaced by H^(x). 
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(4.23) 



(4.24) 




y zy zy z y z 

1a 1b 2a 2b 




y z 
3 



Fig. 1. Feynman diagrams contributing to fx(%o) at one-loop order of 
perturbation theory. The cross indicates the insertion of the axial current 
or density, while the terminal points of the quark lines are at the boundary 
i/o = zo = of the lattice. 

If we now insert the expansion of H(x) in eq. (4.3) and use = 1 and 
(qi_i(x))Q = 0, it follows that 

oo 

fx(x ) = J29o k fxM, (4.26) 

k=0 

/fOro) = ±tT{HW( X yT x HW(x)}, (4.27) 

fi 1 \x ) = ^tv{H^( x yr x H^(x)}) 6 

+ ±( tr{tf( 2 ' ( X yr x H^ (x) + (°) ( X yr x H^ (x)})^ . (4.28) 

The expectation values in the last equation are to be taken at </o = 0. Note 
that (x) is proportional to the gluon field q^y). After contracting the 
gluon fields, the first term on the right hand side of eq. (4.28) may thus be 
represented by the diagram no. 3 in fig. 1. As usual the free quark propagators 
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are depicted by directed solid lines, while the wiggly line stands for the gluon 
propagator. The quark gluon vertex may be inferred from the operator Z)' 1 '. 

The second order correction if' 2 ' (x) is a sum of several contributions. 
Most of them are proportional to q'^ l (y)q h lJ (z) and after contracting these fields 
one ends up with the diagrams no. 1 and 2 in fig. 1. There is one more term, 

HW(x) b = -a^S^H^D^H^iy), (4.29) 

y 

which arises from the 0(a) boundary counterterm (2.4). Since 

SD^^x) = c[ 1] ^{s xo , a [1>{x) - P+^(x - a6)] 

+ 8 X0 , T - a " P-W* + «0)] }, (4.30) 

the integration over the gluon field is trivial in this case. Moreover, using the 
boundary conditions satisfied by the free propagator, it is possible to show that 
the associated contribution, 

/x 1) (^o)b = Re tv{H^{x)^H^ (x) h }, (4.31) 

is of order a (as it should be). 

To compute the Feynman diagrams one inserts the time-momentum rep- 
resentation of the free quark and gluon propagators as given in sect. 3 and 
appendix A. One then has to sum over the time coordinates of the vertices 
and a three-dimensional loop momentum, i.e. there are at most (T/a) 2 (L/a) 3 
terms to be added. Each term involves a few additions and multiplications of 
complex 4x4 matrices. Up to lattice sizes of say 64 X 32 3 , the computational 
work needed to evaluate all diagrams for a given choice of the kinematical pa- 
rameters is then not very large, particularly if the Dirac traces are worked out 
analytically. To guarantee the correctness of our results we have written two 
independent sets of programs (one per author). 
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5. Computation of c A 

From the one-loop results for /a(#o) an d fp( x o) with C = C = it is 
possible to extract the leading term in the expansion 

ca = 4 1) 5o 2 + 4 2) 5o 4 + ... (5.1) 

There are several ways in which this calculation can be done. Essentially 
we shall follow the method outlined in subsect. 6.5 of ref. [1], but to obtain 
a number of cross-checks we decided to approach the problem on a broader 
basis, where one studies the continuum limit and 0(a) improvement of several 
renormalized correlation functions. We shall then also be able to determine the 
one-loop coefficient cj 1 '. 

5.1 Renormalized correlation functions 

To renormalize the correlation functions /a^o) an d fp( x o) we first of all have 
to express the bare parameters too and go through the renormalized ones. We 
adopt the minimal subtraction scheme described in subsect. 3.4 of ref. [1] and 
choose 

to r = 0. (5.2) 
To the order considered the required substitution is then given by 

m =m^gl + O{gi), (5.3) 

g 2 =gi + o(g 4 R ). (5.4) 

For the improved theory the critical bare mass has first been worked out by 
Wohlert [5]. His result, 

N 2 - 1 

oral 11 = -0.2025564(4) x C F , C F = ^ , (5.5) 

has later been reproduced in refs. [8,13] using different methods. 

For the renormalization of /a^o) an d fp( x o) we also need the renormal- 
ization constants Z\, Zp and to one-loop order. So far the latter has only 
been computed in the framework of dimensional regularization [12]. Since the 
renormalization group assumes the same form in all mass-independent renor- 
malization schemes, with coefficients related by finite renormalizations, it is 
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then possible to deduce that 

Z T = -77^FlnM. (5.6) 
The renormalization constants of the axial current and density are given by 

4 1 ' = 0, (5.7) 

4 1} = T^CVlnM. (5.8) 

They have been calculated directly in the lattice theory by considering matrix 
elements between quark states [8,13]. 

The renormalized 0(a) improved amplitudes [/a( 2; o)]r an d [/p( 2; o)]r are 
now defined by 

[/a(zo)]r = Z A Z 2 C {/a(zo) + afsA{xo)} , (5.9) 

[/p(^o)]r = ^p^ c 2 /p(^o). (5.10) 

Note that the factors of the form 1 + &(<7o)ara q need not be written here since 
we are at zero renormalized quark mass which is equivalent to setting m q = 0. 
The correlation function Jsa{ x o) is defined in the same way as /a(#o) with 
A a {x) replaced by the counterterm SA^^x) [cf. eqs. (1.2.26) and (1.6.1)]. To 
one-loop order we have 

\ i „2 / ,(!)/ \ , M) 9 AO), 



+ (z^+2Zf) /fWlfl/UwllOtsU), (5.11) 



[/p (*o)]r = / P 0) (*o) + 5r { /p 1 ' M + m* 1 ) ^-/f 

+ (4 1} +2^) jj 0) ( l0 )} +0(gi), (5.12) 

where all amplitudes on the right hand sides are to be evaluated at too = 0. 
Recall that 

•fx' M = I? (x ) d + A 1] (x ) b (5.13) 
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with fx\ x o)d being the sum of the diagrams no. 1-3 and fx\ x o)b the con- 
tribution of the boundary counterterm proportional to cj 1 '. The explicit 0(a) 
counterterm appearing in eq. (5.11) is given by 

fS M = 4 1 ' \ {dt + d )i 0) (x ). (5.14) 

The renormalized amplitudes are thus numerically known and we may proceed 
to study their behaviour in the continuum limit. 
In the following we set 

T = 2L, O k = 0>O, fi = l/L, (5.15) 

and consider the dimensionless functions 

h P (9,a/L) = [f P (x )] R \ xo=T/2 , (5.16) 

h A (0,a/L) = [f A (x )] R \ xo=T/2 , (5.17) 

h dA (0, a/L) = L\{dt + 9 )[/a(^o)]r| :Co=t/2 . (5.18) 

With properly adjusted coefficients of the 0(a) counterterms all these ampli- 
tudes are expected to converge to the continuum limit with a rate proportional 
to a 2 . To the order of perturbation theory considered, we only require the 
value of c sw at go = which has long been shown to be equal to 1 [4]. The 
coefficients to be determined then are and cj 1 '. 

Taking eqs. (5.11)-(5.14) into account the perturbation expansion of the 
hp , h A and h^A is written in the form 



hp = u + g n \ ui + c| 1) m 2 + am[ 1 K 3 + 



z£ ) +2z£ ) )u \ + 0(gi), (5.19) 



hx = v + g R <, vi + c\ >v 2 + am). >v 3 + 



4 1} + 2Z\} ] ) v + c { l ] Vi + 0(gi), (5.20) 
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hdA = w + gl^w x + c{^w 2 + am^w 3 + 

(4 1} + Wo + ci 1 '^} + 0(g£). (5.21) 

With the exception of u\, v\ and w\, all coefficients appearing in these equa- 
tions can be calculated analytically using the results of sect. 3. In the present 
context we only need to work out their asymptotic form at small a/L. Explicit 
expressions, valid up to corrections of order (a/L) 2 , are listed in appendix B. 

5.2 Computation of to' 1 ' 

Since the renormalized quark mass has been set to zero, the axial current is 
conserved (up to corrections of order a 2 ) and we conclude that 

h dA (6,a/L) = 0(a 2 ). (5.22) 

In particular, the curly bracket in eq. (5.21) is of order a 2 . The same is true 
for the coefficients wq and w 2 , while w 3 grows proportionally to L/a for a — > 
(cf. appendix B). It follows from this that 

am^ ] =-( Wl + c { l ] w 4 )/w 3 + 0(a 3 ). (5.23) 

The term proportional to makes a contribution of order a 2 . The critical 
bare mass may hence be computed by extrapolating the ratio — w\/w 3 to the 
continuum limit. The extrapolation is particularly easy if 9 = 0, because w 4 
is of order a 2 in this case. As a consequence we have 

am^ = -(w 1 /w 3 ) e= o + 0(a 3 ), (5.24) 

i.e. the rate of convergence is improved from order a 2 to order a 3 . 

The only dependence of h^A on the lattice spacing is through the com- 
bination a/L. Taking the continuum limit thus amounts to calculating the 
coefficients Wk for a range of lattice sizes L/a and extrapolating to L/a = oo. 
On general grounds one expects that an asymptotic expansion of the form 

oo 

w 1 /w 3 = a + Y^i a k + b k ln(a/L)](a/L) k (5.25) 

k=l 
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holds [14]. The extrapolation method described in sect. 6 of ref. [14], adapted 
to series with both even and odd powers of a/L, may thus be applied. 

The coefficients Wk (and also Uk and Vk) have been calculated for L/a = 
4, 5, . . . , 32 and three values of 9, equal to 0.0, 0.1 and 1.0. The extrapolation 
of the data at 9 = 0.0 then yields 



This result is consistent with and slightly more accurate than Wohlert's number 
eq. (5.5). Consistent results are also obtained at the other values of 9. 

5.3 Logarithmically divergent terms 

The additive mass renormalization constant determined above cancels the lin- 
early divergent terms in the bare one-loop amplitudes u\, v\ and w\. The 
terms proportional to In (a/L) are cancelled by the renormalization constants 
Z ( l ] , Z { p 1] and Z[ 1] given in subsect. 5.1. From the data for u\ and v\ we can 
now check that the logarithmically divergent terms are indeed as expected. 
From eq. (5.19) we infer that 

Z { p 1] +2Z { ( 1] = constant - ( Ul + u 2 + am { c 1] u 3 ) /u + 0(a 2 ) . (5.27) 

The term proportional to cj 1 ' makes a contribution proportional to a/L to the 
right hand side of this equation. Taking eq. (5.24) into account we obtain 

Zp^ + 2Z^ = constant — (u\ — (w\/ w 3 )g = QU 3 ) / uq + O(a). (5.28) 

It also follows from the above that the error term is reduced to 0(a 2 ) if 9 = 0. 

Before extrapolating the right hand side of eq. (5.28) to the continuum 
limit the constant must be removed by "differentiating" with respect to L, 
i.e. one takes the difference of the elements of the series at L + a and L — a and 
multiplies with —L/2a. In this way one obtains a series which extrapolates 
to the coefficient of the In (a/L) term. The extrapolation then shows that the 
series is compatible with zero as the asymptotic value. A conservative estimate 
of a possible non-zero value is 



am, 



W = -0.2025565(1) X C F . 



(5.26) 



4 1 ' + 2^ 



167T 2 



Cp In (aji) 



c\ < 0.002 



(5.29) 



thus verifying the expected result to high precision. 
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In the same way it is possible to extract Z A + 2Z^ from the one-loop 
contribution to h\. Again it is advantageous to set 9 = and the result 

Z A 1] + 2Z[ 1] = jjr^C F ln(a/x), ~c = -5.9999(1), (5.30) 

is then obtained in complete agreement with eqs. (5. 6), (5. 7). 
5.4 Computation of c A *' 

Since /idAj w o an d w 2 are all of order a 2 it follows that 

c a' = [ _-U7 i + ( w i/ w 3)e=ow 3 ] /w 4 + 0(a), (5.31) 

where 9 should be taken non-zero (-it? 4 vanishes otherwise). The extrapolation 
of the data at 9 = 0.1 then yields 

4 1 ' = -0.00567(1) X C F . (5.32) 

The data at 9 = 1.0 can also be extrapolated with consistent results. It 
should be noted, however, that in this case the order a 2 corrections (which here 
appear as a/L corrections) are almost as large as the leading term. A careful 
extrapolation procedure is then essential to obtain the correct asymptotic value 
of the series. 

It may not be totally obvious that the calculation of presented here 
is in line with the computational strategy discussed in subsect. 6.5 of ref. [1]. 
The proposition made there was to determine ca from the requirement that 
the unrenormalized current quark mass to, defined through eq. (1.6.13), should 
be independent of the kinematical parameters. Since the renormalized quark 
mass has been set to zero, the criterion is that to must vanish up to 0(a 2 ) 
corrections. From the definitions of to and /ija it is straightforward to show 
that 

h d A = Z A Z 2 ( L{2mfp(x ) + ±c A a 3 (9 *9 ) 2 /p M } Xo=t/2 ■ (5-33) 

The last term in this equation is manifestly of order a 3 so that /ija an d m are 
practically proportional to each other. Requiring to to be of order a 2 or /ija 
(as we did here) are hence equivalent conditions. 
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5.5 Computation o/cj 1 ' 

To compute cj 1 ' the best quantity to consider appears to be the ratio 

hp(8,a/L) 

h P (0,a/LY [ ' 

The requirement of 0(a) improvement at one-loop order then implies that 

c| 1 'm2/mo = constant — u\/uq + {ui/uq)b=o 

- om' 11 [u 3 /u - (u 3 /u ) e =o] + 0(a 2 ). (5.35) 

Taking the difference of this equation evaluated at L + a and L — a and multi- 
plying with the appropriate factor then results in a series which extrapolates 
to cj 1 '. Again it is found that the 0(a/L) corrections in this series are not 
small. For both values, 9 = 0.1 and 9 = 1.0, the extrapolation yields 

c[ 1] = -0.0135(1) X C F . (5.36) 

As a further check on this result one may verify that the same number is 
obtained from a similar analysis of the ratio 

h A (9,a/L) 
h A (0,a/L)> 

where the value (5.32) is inserted for . 



6. Non-zero background fields and computation of c s J 

In this section we outline the computation of the first two coefficients in 
the expansion 

c sw = c iJ + c iJ 9o + c sJ 9o + • • • (6-1) 

by demanding the O(a) improvement of the renormalized correlation function 
/idAj defined in subsect. 5.1, for non-zero boundary values C and C . To 
simplify the calculation we drop all O(a) boundary counterterms. This is 
permissible, because the improvement of /ija does not depend on these terms 
(cf. subsect. 6.1 of ref. [1]). 
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6.1 Perturbation theory 

We restrict attention to constant abelian boundary values as specified in sub- 
sect. 6.2 of ref. [1]. For the gauge field V(x,fj,) which minimizes the Wilson 
plaquette action one then obtains [9] 

V(x,0) = l, V(x, k) = exp{aB k (x )}, (6.2) 

B k {x ) = (Cfc - C k )x /T + C k . (6.3) 

This configuration is referred to as the background field in the following. 

The perturbation expansion of h^A is derived essentially as before. The 
main difference is that one expands around the background field, i.e. eq. (4.5) 
is replaced by 

U(x,fj,) = V(x 1 fi) explgoaq^ix)}. (6.4) 

This has already been discussed at length in refs. [9,11] for the case of the pure 
SU(2) gauge theory. Most of the details given there on the gauge fixing and 
the construction of the gluon and ghost propagators and vertices can be taken 
over literally. The extension from SU(2) to SU(3) is straightforward. 

As far as the quark fields are concerned we proceed as in sect. 4. In 
particular, eqs. (4. 2), (4. 3) are still valid and the expansion of ip c \(x) in powers 
of </o is obtained as in subsect. 4.3. For reasons made clear below we now also 
need to expand the quark contribution 

tr {ln(L> + SD + m )} = tr {ln(L> (0) + 5D {0) + m )} 

+ 9o tr{(D^ + 6DW)(dW + 8D^ + mo)" 1 } + . . . (6.5) 

to the total action. The operators occurring in this equation act in the space 
of quark fields with zero boundary values. Note that the first order correction 
does not vanish for non-trivial background fields. 

6.2 Determination of ciw 

Although we already know that ciw = 1 from the original work of Sheik- 
holeslami and Wohlert [4], it is instructive to rederive this result here by ex- 
amining the cutoff dependence of /ija a t tree-level of perturbation theory. The 
calculation can be done analytically in spite of the fact that an explicit expres- 
sion for the quark propagator at </o = is not available. We are only interested 
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in evaluating h<\\ at zero quark mass and thus set too = throughout this sub- 
section. 

As in the case of vanishing boundary values we have 

/f (x ) = -±tT{HW(xy 7o Hl°\x)}, (6.6) 

where the matrix if' ' (x) is defined through 

(£)(o) + 8D^)H^ (x) = 0, < x < T, (6.7) 

P + hW(x)\ X0=0 = P+, P-HW(x)\ Xo=t = 0. (6.8) 

(cf. sects. 2 and 4). The solution of these equations has a relatively simple 
form as a result of the special properties of the background field. First note 
that 

B k (x ) = ib(x ) (6.9) 

for all k = 1,2,3, where b(xo) is a linear function of xq with coefficients that are 
real diagonal matrices in colour space. In particular, the background field is 
invariant under space translations and the solution if' ' (x) is hence a function 
of the time xq only. Eq. (6.7) thus becomes 

{P + dt - P-d + A(x ) + iB(x )j + iC(x ) lol }H^(x) = 0, (6.10) 

where j = ji + 72 + 73 and 

A(x ) = ^ {sin [\a{b{x ) + e/L)] f , (6.11) 
B(x ) = ^-sin[a(b(x ) + 9/L)], (6.12) 

(0) 

C(x ) = -^-sm[a 2 d b(x )]. (6.13) 
To solve eq. (6.10) we make the ansatz 

HW(x)= [s(x ) + it(x )j]P+. (6.14) 
The Dirac equation and the boundary conditions then reduce to 

dos = 3(B + C)t- As, s(0) = 1, (6.15) 
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d t = (B-C)s + At, t(T) = 0. 



(6.16) 



Note that the coefficients A, B, C and the amplitudes s and t are all diagonal 
matrices in colour space with real entries. The solution of the equations can 
be worked out numerically, but this will not be required in the following. 
From eqs. (6.6) and (6.14) we now deduce that 

ff ) (x ) = tr{-s(x ) 2 + 3t(x ) 2 }, (6.17) 

where the trace is over the colour indices. Using the equations above it is then 
straightforward to show that 

1 1 (d* + d )f { A 0) (x ) = 

tr{ci(a;o)s(a;o) 2 + c 2 (x )t(x ) 2 + c 3 (xo)s(x )t(x )} (6.18) 

with coefficients c\, c 2 and c 3 given explicitly in terms of A, B and C. In 
particular, their behaviour in the continuum limit is analytically calculable 
and one finds that c\ and c 2 are of order a 2 while 

c 3 (x ) = a(4°J - l)d b(x ) + 0(a 2 ). (6.19) 

The amplitudes s and t converge to some well-defined smooth functions in the 
continuum limit and we thus conclude that 

4°i = a (4°J - l)Lti{dob(xo)s(xo)t(x )} Xo=T/2 +0(a 2 ). (6.20) 

The trace can be worked out explicitly using the formulae given in appendix 
C and it is then immediate that it does not vanish for general boundary values 
of the specified type. For 0(a) improvement the coefficient c'w must hence be 
set to 1. 

6.3 Computation of c sw to one-loop order 

The one-loop diagrams contributing to fx( x o) are listed in figs. 1 and 2. The 
lines in these diagrams now represent the gluon, quark and ghost propagators in 
the given background field. The tadpole diagrams in fig. 2 with the quark loop 
arise from the combination of the first order term in the expansion (6.5) with 
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y z 

4a 



y z 

4b 



= -N^N^ \ -|- 



Fig. 2. Tadpole diagrams contributing to fx(xo)- The directed dashed 
line represents the ghost propagator. 

the first order correction to tv{H (x)^TxH (x)} . These are the only diagrams 
proportional to the number N{ of quarks. 

Since there is no general closed analytic expression for the propagators in 
the background held (6. 2), (6. 3), the computer time required for the calculation 
of the diagrams is not negligible. We have hence decided to set T = L (rather 
than T = 2L) and to produce data only for even L/a up to L/a = 32. As for 
the boundary values of the gauge held we made the same choices as in ref. [10]. 
We now also set c'w =1, = or = 1 and perform all calculations at zero 
quark mass. 

The perturbation expansion of h^A to order g R assumes the form 
hdA = w + #r jwi + am^w 3 + (^Z^ + w 

+ c { l ) w 4 + 4 1 Jw 5 ^+0(g 4 R ). (6.21) 

The coefficients wq, w\, . . . are not the same as in eq. (5.21), but their definitions 
are completely analogous. Since wq is of order a 2 the corresponding term in 
the curly bracket may be dropped in the following. From the discussion in 
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subsect. 6.2 one deduces that 



w 5 = aLtr{dob(x )s(xo)t(x )} xo=T/2 + 0(a 2 ). (6.22) 

For the chosen boundary values of the gauge field, the trace in this expression 
does not vanish. We thus conclude that 

4$ = "K + am^w 3 + c { l ] w 4 )/w 5 + 0(a/L). (6.23) 

the valu 

yields 



Using the value of given in eq. (5.32), the extrapolation of our data then 



0.155(1) for N = 2, 
4J = { (6-24) 
\ 0.267(1) for TV = 3. 

There is no dependence on N{, because the contribution to w\ of the quark 
tadpole diagrams turns out to be of order a 2 . 

Our computations have been checked in various ways. For N = 2 the 
expressions for the ghost and gluon tadpoles can be compared with the cor- 
responding expressions given in ref. [11]. We have also verified that the total 
one-loop contributions to /a(#o) an d fp( x o) are independent of the gauge pa- 
rameter. For this the gluon tadpole diagrams are essential. In addition runs 
up to L/a = 20 were made with alternative choices of the boundary values C 
and C . In all cases the values quoted above for ciw have been reproduced, 
although to a reduced precision. Finally by combining our data at 9 = and 
= 1, it is possible to extract the coefficient . The values that have been 
obtained in this way are less accurate but completely consistent with eq. (5.32). 

Our results confirm those obtained by Wohlert [5] nearly 10 years ago, 
who actually cites numbers with a higher accuracy, viz. 

f 0.15400(4) for N = 2, 

(6.25) 

( 0.26590(7) for TV = 3. 

We are unfortunately unable to directly check whether Wohlert 's error analysis 
is realistic since his raw data are no longer available. But our result, eq.(6.24), 
is perfectly satisfactory for application at small bare coupling. At higher values 
of </o one should in any case compute c sw non-perturbatively [2]. 
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Table 1. Values of the coefficients ro and r\ at mo = m c , T = 2L and xo = T/2 



a 

U 


T.ln 

1 J 1 It 


1 


1 1 


' 1 \c W -0 
o A — u 


0.0 


6 


0.0 


0.000861 


0.000861 




8 


0.0 


0.000175 


0.000175 




10 


0.0 


0.000055 


0.000055 




12 


0.0 


0.000023 


0.000023 




14 


0.0 


0.000011 


0.000011 




16 


0.0 


0.000006 


0.000006 


1.0 


6 


0.001151 


0.001027 


0.001963 




8 


0.000365 


0.000254 


0.000782 




10 


0.000150 


0.000092 


0.000431 




12 


0.000072 


0.000042 


0.000278 




14 


0.000039 


0.000022 


0.000195 




16 


0.000023 


0.000013 


0.000146 



7. How large are the remaining cutoff effects? 

To answer this question we here consider a few quantities constructed from 
/a an d fp an d study their approach to the continuum limit. For simplicity we 
restrict attention to the case of vanishing boundary values of the gauge field 
(the situation is not very different at non-zero values of C and C"). 

7.1 PC AC relation 

As in sect. 5 we again set T = 2L, 9 k = 9 and study the theory at the critical 
point mo = m c . The unrenormalized current quark mass to, defined through 
eq. (1.6.13), is then expected to vanish up to terms of order a 2 . The coefficients 
ro and r\ in the expansion 

am = r + riCp^o + O(5o) (7.1) 

should hence go to zero in the continuum limit with a rate roughly proportional 
to (a/L) 3 . The numbers listed in table 1 show that at xq = T/2 both, ro and 
ri, are indeed small and rapidly approaching zero. 
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It is interesting to note in this connection that the contribution of the 0(a) 
correction proportional to is negligible at small values of 9, while at 9 = 1 
a visible improvement is achieved by including this term. The contribution of 
the boundary counterterm proportional to cj 1 ', on the other hand, appears to 
be zero within rounding errors and is certainly completely negligible. 

The values of am that one obtains from eq. (7.1) and the coefficients listed 
in table 1 are quite small. At L/a = 8, 9 = 1, g$ = 1 and N = 3, for example, 
the value am = 0.0007 is found. In the quenched approximation, where the 
cutoff l/a is about 2 GeV at g$ = 1, this translates to a mass m ~ 1.4 MeV. 
As discussed in subsect. 6.6 of ref. [1], this may be taken as an indication for 
the systematic uncertainty which one has when defining the zero mass point 
on a 16 X 8 3 lattice. By increasing the lattice size to say L/a = 12 (which is 
still rather small by normal standards) the uncertainty is reduced by a factor 
6 at 9 = 1 and nearly an order of magnitude at 9 = 0. One should however 
be warned that perturbation theory can only give a rough idea on the actual 
size of the cutoff effects. In particular, at the larger couplings the lattice effects 
associated with the non-perturbative scales may become important (cf. ref. [2]). 

7.2 Another test of PCAC 

So far we have set xq = T/2 in the definition of m, but up to 0(a 2 ) corrections 
the same value of m should be obtained for all xq. As discussed in sect. 3 the 
lowest order coefficient ro is in fact independent of xq. 

The behaviour of r\ as a function of xq on a 16 X 8 3 lattice is shown in 
fig. 3. A significant increase of the cutoff effects towards the boundaries of 
the lattice is observed, but it should also be said that they remain small and 
certainly are not above the expected order of magnitude. The fact that the 
effects are larger near the boundaries may have several causes. An obvious 
possibility is that the contributions of the high-energy states to the correlation 
functions /a^o) an d fp( x o) are n °t sufficiently suppressed if xq or T — xq is 
small in lattice units. But it could also be that the effects seen in fig. 3 arise 
from the 0(a 2 ) corrections associated with the low-energy intermediate states. 

On larger lattices the cutoff effects are rapidly becoming small and when 
r\ is plotted against xq the picture is essentially a scaled version of fig. 3. An 
interesting detail is that r\ is practically independent of 9 in the middle of 
the lattice. For this to come out it has been important to include the O(a) 
correction proportional to , which is independent of xo and approximately 
given by 6c^' (a9/L) 2 . 
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Fig. 3. Values of r\ on a 16 x 8 3 lattice as a function of the time xo at 
which the axial current is inserted. Circles and crosses are for = and 
= 1 respectively. 

7.3 Cutoff dependence of fp(xo) 

The definition of the renormalized correlation function [/p(:eo)]r. involves some 
arbitrariness in the choice of normalization. It may then not be easy to decide 
whether the cutoff effects that one finds should be considered large or small. 
More sensible quantities to study are ratios of correlation functions where the 
normalization factors cancel. 

To illustrate this we introduce the function 



F(9, a/L, g ) = f P (x ) \ Xo=T /2 , 
where T = 2L, mo = m c and 9k = as before. The ratio 

F(9,a/L',g ) 



P 



F(0,a/L,g o ) 



L' = 2L, 



(7.2) 



(7.3) 



then is a renormalization group invariant, i.e. in the continuum limit p con- 
verges to a function of 9 and the running coupling g 2 at scale L (in any scheme). 
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Table 2. Values of the coefficients po and p\ at mo = m c , T = 2L and «o = T/2 





T In 


Po 


Pi 


Pi\- C w =0 


0.0 


6 


1.0 


0.0236 


0.0236 




8 


1.0 


0.0101 


0.0101 




10 


1.0 


0.0056 


0.0056 




12 


1.0 


0.0036 


0.0036 




14 


1.0 


0.0025 


0.0025 




16 


1.0 


0.0019 


0.0019 


1.0 


6 


0.9584 


-0.0049 


-0.0099 




8 


0.9755 


-0.0033 


-0.0076 




10 


0.9839 


-0.0022 


-0.0059 




12 


0.9886 


-0.0016 


-0.0048 




14 


0.9916 


-0.0011 


-0.0040 




16 


0.9935 


-0.0008 


-0.0034 



The approach to the continuum limit of the ratio (7.3) can now be studied 
in perturbation theory. In table 2 the coefficients po and p\ appearing in the 
expansion 

P = Po + PiC F g 2 + 0(g 4 ) (7.4) 

are listed for two values of 9 and a range of L/a. In the continuum theory the 
coefficients are 1 and 0, respectively, for all 9. The cutoff effects thus appear 
to be rather small and decrease roughly proportionally to (a/L) 2 as expected. 
From the third column in the table we can also see that the O(a) counterterm 
proportional to cj 1 ' is completely negligible at 9 = 0. When 9 increases the 
contribution of the boundary counterterm can become numerically significant 
compared to the cutoff effects at tree-level, but it remains rather small in 
absolute terms if (say) Cpg 2 < 1. 
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Appendix A 

The Fourier representation of the gluon propagator (4.17) reads 

D^(x, y) = L- 3 J2 eip{x ~ y)d ^(xo, yo; P). (A.l) 
p 

For any momentum p we define the "energy" e > through 

cosh(ae) = 1 + |a 2 p 2 . (A. 2) 

In view of the symmetry 

d^{x , y ; p) = d UfJi (y , x ; p), (A. 3) 

it suffices to consider the case xq > yo- In the Feynman gauge, Ao = 1, and for 
non-zero momenta p the propagator is then given by 

doo{x , y ; p) 



sinh(ea) sinh(eT) 
X cosh[e(T — x — \a)] cosh[e(y + \o)'\, (A-4) 



dkj (xo , yo ; p) = hj -T-r-7 — s . i (rr , 
smhfea) smhfei J 

X sinh[e(T - x )] sinh(ey ), (A. 5) 

while for p = we have 

doo{x , y ; 0) = y + a, (A. 6) 

dkj(x , yo] 0) = S kj (T - xq)^. (A. 7) 

The mixed components c?ofc an d ^fco vanish for all time coordinates and all 
momenta p. 
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Appendix B 



Most of the coefficients appearing in eqs. (5.1 9)-(5.21 ) are known analyti- 
cally and can be expanded in powers of a/L. It is convenient to introduce the 
abbreviations 

si = sinh (2y/30), (B.l) 
co = cosn (2y/30). (B.2) 

Up to corrections of order (a/L) 2 the analytically calculable coefficients are 
then given by 

uo = — , (B.3) 

CO 

127V6»si a . 
Nsi L ( 4N0 2 (l - si 2 ) N0si 1 a 



U3 V30co 2 a + \ 3co 3 ^co 2 )^ (B ' 5) 



= — ^ , (B.6) 

2AN0s\ a 

^ = 7- B.7 



\/3co 3 L 
TVsifco-2) L f JV0si(co-2) 



V 3 



^/30co 3 a { 6^3 co 3 

AN0 2 , o 9 , 1 a , , 

+ ^cT ( ° + 6) jL' ( 8) 

67V#si a _ . 



w = 0, (B.10) 
w 2 = 0, (B.ll) 
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2N L f 8N0 3 si 3N0 2 ) a , 

^3 = + -j=— + 7 , B.12 

co a y a/3co 2 co J L 

12N0 2 a 

w 4 = -. B.13 

co L 



Appendix C 

The functions s and t introduced in subsect. 6.2 can be computed analyt- 
ically in the continuum limit. The equations to be solved are 

d s = 2(b + 6/L)t, s(0) = 1, (C.l) 

d t = (b + 9/L)s, t(T) = 0. (C.2) 

To be able to write down the solution in a compact form it is useful to define 

Note that the denominator dob(xo) is independent of xq. One may then easily 
verify that 

s(xq) = ky/s cosh | [^(xo) — z(T)] , (C4) 
t(xo) = k sinh | [^(a;o) — z(T)] , (C5) 
solves the system (C.l), (C.2) if we choose 

k = |V3coshi[^(0) - z(T)]} 1 . (C.6) 
It is also not difficult to show that the solution is unique. 
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